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' Abstract. In this paper we study a generalization of the ID Hubbard model by considering 

' density-density and Ising-type spin-spin nearest neighbor (NN) interactions, parameterized by 

V and J, respectively. We present the T = Q phase diagram for both ferro (J > 0) and 
' anti-ferro (J < 0) coupling obtained in the narrow-band limit by means of an extension to zero- 

temperature of the transfer-matrix method. Based on the values of the Hamiltonian parameters, 
we identify a number of phases that involve orderings of the double occupancy, NN density and 
' O ' spin correlations, being these latter very fragile. 
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1. Introduction 

The Hubbard model (HM) is undoubtedly one of the most studied model in condensed matter 
\ physics. In one spatial dimension (ID), the exact Bethe solution for the pure HM is known [1]. 

CNl ■ During its long history, several extensions to the HM have been proposed. Among them, the 

inclusion of non-local charge density and/or spin exchange interactions appears to be the most 
natural choice. Density correlations correspond to an effective finite-range Coulomb interaction 
(repulsive or attractive) as opposed to the on-site Hubbard U . On the other hand, the spin-spin 
interaction is believed to be responsible for the interplay between various spin orderings (SDW 
or ferromagnetism) with strong (triplet) superconducting correlations in a variety of compounds: 
. Sr2Ru04 [2], UGea |3|, URhGe @], ZrZng [5] and (TMTSF)2X [6] family. Any extension of 

r> \ the HM breaks the integrability, so that either approximate analytic {e.g. bosonization) or 

' numerical finite-system methods {e.g. Density Matrix Renormalization Group or Quantum Monte 

Carlo) should be used even in ID. Bosonization works well in the limit of large bandwidth 
and predicts several interesting effects like coexistence of spin density wave correlations and 
(triplet) superconducting ones (for spin exchange interactions [3 |5]) or a transition towards 
Mott-insulating phases at commensurate fillings (for charge interactions [9]). On the other 
hand, in the opposite limit of narrow-band width (atomic limit), the extended Hubbard model 
can be exactly solved by means of any of the two equivalent methods - Transfer Matrix (TM) 
technique |10| or Composite Operator Method |1H [T2] . In the present work we employ a T ^> 
extrapolation of the TM method in order to derive the zero-temperature phase diagram of the 
extended HM with both nearest-neighbor (NN) charge and spin interactions in the narrow-band 
limit. 



2. Results 

The Hamiltonian of the U-V-J system reads as foUows: 



H = Y^ [-/^"W + UD{i)] + y ^ n(i)n(i + 1) - J J]] 5^(i)S^(i + 1). (1) 

i i i 

Here n{i) = + n^(i) denotes the electron density operator at site D[i) = n^(i)n|(i) is the 
double occupancy at site i while S^{i) = {n^{i) — nj,(i))/2 denotes the 2-component of the spin 
at site i. The Hamiltonian ([TJ contains only NN term, hence the TM will be a 4 x 4 matrix. 
It is easy to check that the TM elements can be calculated through the matrix elements of an 
auxiliary matrix Z: Tij = exp (—fiZij), where Z is defined as follows: 



Z 



1 V-{-fi V + ^-ii I +'21/- I/. 



(2) 



In our longer paper |13| we show how to reconstruct the T = phase diagram from the TM 
matrix elements. Here we apply this method to the U-V-J model. We first identify all different 
matrix elements (we call them energy scales). In the present case these are: 

=0 

Fl =u^-^,- F^ = V-{-^l■ Fl = V + {-^^ (3) 
F3/2 =ll + 2V-y 
=u + W - 2fi. 

In ([3]) we have already sorted the energy scales based on the values of particle numbers (given 
by the superscript of F) allowed for a given configuration. In particular, the three energy terms 
at n = 1 can be written as a unique expression as follows: 

F' = A-i,, (4) 

where: 

/U \J\\ { IL ^ T/ _ M 

V2' 4; f>^-^. ^' 

We choose J as energy scale and explore the phase diagram in the V—U plane. As fi increases, the 
particle number in the system increases as well. Therefore, for each value of n (n = 0, ^, 1, |, 2) we 
can establish the ranges of the chemical potential, within which fi can change without changing 
n. Such ranges can be summarized as follows: 

n = : fl < min (O, A, x, |x) n = 1 : max {A, 2A) < fi < 2x — max {A, 2A) 

n = \ : < /i < min {2A, x, ^x) n = | : max x, 2(x — Aj) < < 2x (6) 

n = 2 : max (x, |x, 2x — A, 2x) < fi, 

where we introduced x = ^(C/ + 4V). Except for the cases n = and n = 2, which are always 
possible, the other cases might be incompatible with a given choice of U, V, J. This occurs 
when the upper-bound of ^ becomes lower than the lower-bound. This condition determines the 
existence of the jump at a given n in the dependence fi{n). Such conditions are: 

n=i X > A A> 

n=l0<2A<xoi:A<0AA<x (7) 
n=| X > A A> 0. 
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Figure 1. Left panel: Phase diagram in the case < n < ^. Right paneh Phase diagram in the 
case 4 < n < 1. 



The conditions for the existence of jumps at n = ^ and n = | coincide as these values of n 
are related by the particle-hole relation. Moreover, depending on the values of the Hamiltonian 
parameters, we can distinguish four cases in the conditions ([6|): 

(i) x > A ^ > (ii) x > A ^ < (iii) x < A yl > x (iv) x < A ^ < x. 

We can identify three zones in the V — U plane depending on how many jumps there are in 
the dependence //(n): 

(i) y > ^ A [/ > 

(ii) (a)f>-1/-^Af>F-^Ay<^ (b)C/<OAT/>OAf<F-^ 

(iii) F < A [/ < -21/- ^. 

The zone i) is characterized by the presence of jumps at all "allowed" values of n: n = 0, ^, 1, |, 2, 
while in the zone ii) the jumps are only at n = 0, 1, 2. Finally, in the zone iii), the only jumps 
present are at n = 0, 2. The dependence /-i(?i) in the three zones can be summarized as follows: 



n = 
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iii) 

We are now in a position to transform from canonical (n fixed) to grand canonical (^ fixed) 
ensemble. In the zone i) and in the range 0<n< ^,/x = and the free-energy per site 
(which is a linear function of n) goes from Fi = to -F2 = — /u/2|^=o = 0. Therefore, in the 
above interval F is constantly zero and hence the internal energy per site, which is defined as 
E = [{H) + fJ-^i{n{i))] = F + fin, is zero as well. Analogously, we can determine the behavior 
of the internal energy as a function of n in the whole range. These findings for the zones i) and 
ii) are summarized in Tables [T] and [21 In the zone iii), it is easy to verify that E = (U + 4:V)^ 
for all < n < 2. 



Table 1. Internal energy as a function of particle concentration E{n) for the Zone i). 



< n < i 


5 < n < 1 


1 < n < 1 


f < n < 2 


E = 


E = A{2n - 1) 


E = {U + W){n - 1) + A{3 - 2n) 


E = {U + 4V){n - 1) 



Table 2. Internal energy as a function of particle concentration E{n) for the Zone ii). 



< ra < 1 


1 <n<2 


E = An 


E={U + 4V){n - 1) + A{2 - n) 



By definition, A takes two different forms depending on whether ^ is greater or not than 
y - ^. This means that everywhere A appears, the two cases should be considered and, 

therefore, the zones i) and ii) will be further divided in two by the line ^ = V — Thanks to 
the particle-hole symmetry enjoyed by the system, we can consider the phase diagrams just for 
n ^ 1. Two cases can be distinguished: < n < ^ and ^ < n < 1, together with the particular 
cases n = ^ and n = 1. The two phase diagrams are depicted in Fig. [1] We proceed now to the 
description of the different phases appearing in the whole phase diagram. 

Phase a is located in the range V > AU > and < n < ^. The internal energy is zero 
because the particles are are far apart to activate any of the Hamiltonian terms. The chemical 
potential vanishes, /i = 0. When n > ^, the former Phase a gives rise to two new phases (?/ and 

C). 

Phase (3 is located in the range ^ > \V\ — A F < Energy: E = {V — and thus 

in this phase the singly occupied sites at NN distance interact via density-density and spin-spin 
interactions. The spin-spin correlation function takes the value: 

{5W(i+l)> = -§f = 24^. (8) 

Hence, if J < 0, the system is dominated by antiferromagnetic correlations, while if J > the 
correlations are ferromagnetic. This phase exists in the range < ?i < 1. Chemical potential 
;U = y — NN charge-density correlations {n(i)n(i + l)) = n. This phase exhibits spontaneous 
magnetization, uniform in case J > 0, or staggered in case J < 0. 

Phase 7 is located in the range V > A U <OA-^<y — Energy is E = ^n, and 
therefore in this phase the doubly occupied sites are separated by at least one empty site. This 
phase exists in the range < n < 1. Chemical potential fj, = ^. The only non-zero correlation 
function is {D{i)) = n/2. 

Phase 6 is located in the range V < OA^ < —V — for all < n < 2. The energy: 
E = (U + 4V)^ and thus, in this phase, the doubly occupied sites are placed at the NN distance. 
Chemical potential fj, = ^ -\- 2V. The only non-zero correlation functions are {D[i)) = n/2, 
{n{i)n{i + 1)) = 2n and {D{i)D{i + 1)) = n/2. 

Phase rj is located in the range V > -^A^ > V — and ^ < n < 1. The 

energy E = (V — ^)(2n — 1). This is the second phase with magnetic correlations and is 
similar to the adjacent Phase f3 in the sense that in both phases density-density and spin-spin 
correlations are present, although the expressions for the correlation functions are different: 
{n{i)n{i + 1)) = 2n - 1, (5^(i)S'^(i + 1)) = 2iS^(2n - 1). Chemical potential fi = 2{V - ^). 
There also exists spontaneous magnetization, uniform in case J > 0, or staggered in case J < 0. 

Phase is located in the range y > V —^-^AU > and ^ < n < 1. The energy E = U{n—^). 
This phase is similar to the Phase 7 in the sense explained above. Chemical potential ^ = U, 
while the only non- vanishing correlation function is {D{i)) = n — 1/2. 



The phase boundaries for the phases at n > 1 are identical to those at n < 1 and the 
expressions for the internal energy can be obtained from those at n < 1 by invoking the particle- 
hole symmetry: 

E{2-n) = E{n) + {U + W){l-n), (9) 
while for the chemical potential the particle-hole relation states: 

fi{2-n) = {U + W) - fi{n). (10) 



3. Conclusions 

In this work we present an original method for extrapolating the TM technique at T = and use 
it to obtain the exact T = phase diagram of the extended Hubbard U-V-J model in the narrow- 
band limit. Depending on the values of the Hamiltonian parameters, the orderings involve double 
occupancy, NN density and spin correlations. We note that the phase diagram boundaries do 
not depend on the sign of J, although the spin correlations do. The spin correlations appear 
to be rather fragile in our model: the Ising term in the Hamiltonian acts only if the NN sites 
are singly occupied. However, in the same conditions the V term is also active, while the U 
term tends to create either doubly occupied or empty sites (depending on the sign of U). That 
is why, when either V is large and positive or U is large and negative the spin correlations are 
completely suppressed. The introduction of J affects only the phases (/3 and rj) where the sites 
are singly occupied. In these phases, for J = |14| I12|. the spins are not ordered while, if J 7^ 0, 
a ferromagnetic or antiferromagnetic order is established depending on the sign of J. 
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